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Abstract 

The formation of singularities in the three-dimensional Euler equation is 

investigated. This is done by restricting the number of Fourier modes to a set 

which allows only for local interactions in wave number space. Starting from 

an initial large-scale energy distribution, the energy rushes towards smaller 

scales, forming a universal front independent of initial conditions. The front 

results in a singularity of the vorticity in finite time, and has scaling form 

as function of the time difference from the singularity. Using a simplified 

model, we compute the values of the exponents and the shape of the front 

analytically. The results are in good agreement with numerical simulations. 
PACS-numbers: 47.15Ki, 47.27Eq, 47.54-Fr 
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I. INTRODUCTION 



The aim of the theory of fully developed turbulence is to understand fluid flow at very 
high Reynolds numbers. Energy which is fed into the system at some outer scale L is 
transported to increasingly smaller scales through a series of instabilities, until this cascade 
is stopped by the smoothing effect of viscosity. It therefore seems natural to consider the 
limiting equation where the viscosity is put to zero, and the Reynolds number thus infinite. 
The resulting Euler equation will not be able to describe a stationary state, where the 
influx of energy is balanced by viscous dissipation. Rather, the expectation is that as 
viscosity no longer limits the smallest excitable scale, the breakdown of structures will 
continue indefinitely and a singularity of the derivatives of the velocity field will appear in 
finite time, as first suggested by Onsager 

This singularity has attracted considerable attention, in particular from a numerical 
point of view . The reason is twofold: First, the singularity is a refiection 

of the instability of turbulent structures, and thus should give insight into the mechanism 
for energy transfer in fully developed turbulent fiow. Second, as the spatial and temporal 
scale of the singular fiow gets smaller, one expects the solution to become independent of 
boundary or initial conditions, and to refiect the properties of the nonlinear interaction 
alone |]IU|. The resulting solution or class of solutions should thus represent a "coherent 
structure" of the turbulent fiow as well, as long as viscosity is not yet important on the scale 
of its spatial variation. Such a structure is a likely candidate to represent the small scale 
structure of a turbulent velocity field, which has become independent of its outer boundary 
conditions. 

But despite considerable numerical efforts, singularities of the Euler equation have re- 
mained elusive. There is disagreement about their expected structure, and even the very 
existence of a singularity is a subject of debate. Previous papers have been about equally 
divided between giving indications in favor p|,||,|^J^ or against the existence of a 

singularity. In this paper we propose to attack the problem using cascade models, which 
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have been widely used to study fully developed turbulence. The idea of cascade or shell 
models is to divide wavenumber space into bands, which cover a certain ratio in wavenum- 
ber. Between different bands only local interactions are permitted, thus implementing the 
physical idea of local transfer originally proposed by Kolmogorov [|1T|. This results in a 
tremendous simplification of the problem, both conceptually and numerically. We will thus 
be able to confirm the existence of a singularity unambiguously and to study its properties 
in great detail. Moreover, further simplification of the model will allow us to find analytical 
solutions, which confirm the existence of a unique singular shape. This is particularly useful 
since it provides a unified description of both Euler and Navier-Stokes dynamics in terms 
of a single cascade model. The solid understanding of both aspects should enable us to 
ascertain the significance of Euler singularities to turbulent flow. 

In the next section we will introduce a class of models originally developed in known 
as REduced Wave vector set Approximations or REWA models. They arise by restricting 
the number of available Fourier modes to a self-similar set, with a constant number of 
modes within an octave in wave-number. The properties of REWA models have been studied 
extensively in the context of stationary turbulent flow |T2|,|TB|,|T^ . In particular, it was shown 



that the turbulent fluctuations are characterized by a set of anomalous scaling exponents. 



as suggested by the multifractal theory of turbulence [jT5|. Here we make the connection 



between the inviscid singularity and the stationary state of turbulent flow by presenting a 
simulation of decaying turbulence, starting from an initial large-scale distribution of energy. 
If the viscosity is sufficiently small, the flow will be effectively inviscid, resulting in a rapid 
build-up of velocity gradients. Eventually, after sufficiently small scales are excited, viscosity 
becomes important, leading to dissipation of energy. As inertial transport and viscous 
damping balance, the energy spectrum becomes flatter and close to a Kolmogorov spectrum. 

In the third section we study the formation of singularities for very long cascades at 
zero viscosity. Starting from arbitrary initial conditions, a universal front develops, which is 
self-similar: at different time distances from the singularity the solutions can be collapsed 
by a rescaling of their length scale. The smallest excited scale Vc follows a power law as 



function of the distance t = t* — t from the singularity: Vc ~ r^. The relevant exponents 
and the form of the front are determined. In the fourth section we develop an analytical 
description of the singularity by using an effective equation for the energy of the shell. The 



same effective equation has been used before |T6| , p!^ to compute the anomalous exponents 
of stationary turbulence. We show that fluctuations are irrelevant for the description of 
the Euler singularity. The resulting deterministic equation can be reduced to an ordinary 
differential equation if the self-similarity of singular solutions is exploited. 

This ordinary differential equation is used in the fifth section to study the selection of 
universal solutions out of arbitrary initial data. There exists a family of solutions parame- 
terized by the exponent /?, which connects length and time scales. For large (3, the solutions 
develop unstable fronts, which are unphysical. Thus the most singular solution which is 
not yet unstable is selected. The resulting unique solution agrees well with numerical sim- 
ulations of the REWA cascade. In the discussion we comment on related work and point 
out possible uses of the present study of inviscid singularities for the understanding of fully 
developed turbulence. 



II. MODEL EQUATIONS 

A variety of shell models have been proposed in the past to describe a turbulent cascade 
||T7| , p^JT9| , |20|Jl^ . By using a dynamical model, one hopes to gain insight into the origin and 
the statistics of turbulent fluctuations. Since a cascade model consists of a linear structure 
of turbulence elements, the problem is simplified enormously, both from an analytical and a 
computational point of view. However, investigations of cascade models have been limited 
to the steady state, where energy input equals dissipation on the average. The question ad- 
dressed here is whether cascade models are also capable of describing some of the instabilities 
of inviscid flow, which lead to the build-up of gradients. 

The cascade model we consider here was introduced in [jl2|, and is sometimes called the 
REWA model. It has been used extensively to study the stationary state of fully developed 
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turbulence ||T^,|T3|JT^ . It is based on the full Fourier-transformed Navier-Stokes equation 
with a volume of periodicity {27tL)^. Only local interactions are taken into account, which 
is implemented by projecting the Navier-Stokes equation onto a self-similar set of wave 
numbers /C = f(^e- Each of the wave vector shells JCe represents an octave in wave number, 
which greatly reduces the total number of modes, making the model numerically tractable. 
The shell /Co describes the motion of the largest elements in the flow, which are of the order 
of the outer length L. It is composed of wave vectors kf ^ : /Co = {kf ^ : i = 1, . . . , n}. 
Starting with the generating shell /Co, the other shells are found rescaling /Co with a factor 
of 2: ICe = 2^/Co. The shell JCe thus represents structures of size r ~ 2~^L. In a turbulent 
cascade, this scaling procedure is followed until one reaches a Kolmogorov length rj, where 
the turbulent motion is damped by viscosity. In the present paper, we will mostly be 
concerned with the limit of zero viscosity. Thus arbitrarily small scales can be excited, and 
our simulations are valid only for a finite time, until energy is transferred into the smallest 
scale available. By choosing the number of levels very large, we are still able to extract 
reliable scaling information. 

Explicitly, the projection of the Navier-Stokes equation reads 
d 

—Ui{k,t) = -iMijkik) Uj{p,t)uk{<i,t) - uPui(k,t) (la) 

k=p+q 

k-u(k,t) = . (lb) 

The coupling tensor Mjjfc(k) = [kjPik{k) + kkPijik)] /2 with the projector Pjfc(k) = 6ik — 
kikk/k"^ is symmetric in j,k. The inertial part of ([Ta| ) consists of all triadic interactions 
modes with k = p + q. For the moment we have kept the viscous term, but will put z/ = 
later for our study of the Euler equation. With this approximation the energy of a shell is 

Ei{t) = li: |u(k,t)r . (2) 
Before the appearance of a singularity the total energy 

oo 

Etotit) = J2Eeit) (3) 

£=0 



is conserved for z/ = 0. 



In an earlier paper we have investigated the effect of different choices for the wave 



vector set /Co in some detail Here we will mostly deal with a single set of = 26 modes, 
where the components of k-'^'* consist of all combinations of 0, —1, and 1, because we found 
the properties of the inviscid singularity to be quite insensitive to the specific choice of /Cq- 
The /Co considered here only allows for local couplings between shells, and thus the energy 
is transported only between adjacent shells. Hence if T^^^+i(t) is the energy transfer from 
shell i to i + 1, and no dissipation occurs, we can write an energy balance equation 

j^Ee{t) = Te_,^eit) - Te^e+iit) . (4) 

The transfer T^_,^+i(i(;) can be written explicitly as a sum over triple products of velocity 
modes. Equation (|^) will be the basis for a simplified description of the cascade, which we 
will use later to obtain analytical solutions. 

We illustrate the formation of a singularity in our model by considering the dynamics 
(|I|) with an initial condition where only the modes of level i = are excited. To make 
a connection with stationary turbulence, we keep u small but finite in this example. The 
Reynolds number 

LU 

Re = , 

u 

where U is the typical amplitude of a velocity mode on the highest level, is 1.35- 10^. Figure 
ID shows the resulting evolution of the shell energies. 
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FIG. 1. Evolution of the REWA cascade at a Reynolds number of 1.35 • 10®. The energy is 
initially localized in the 1 = level. Viscous effects are small at first, and a scaling regime with 
exponent a = 1.12 develops. After the singular front is stopped by viscosity, the energy levels off 
to a Kolmogorov spectrum. 

Energy rushes downward in scale to fill the shells which are not yet excited. As shown in 
TM, this can be seen as a result of the tendency of the dynamics to establish equipartition 



of energy between shells. Since the small scales are not excited at all, the energy transfer 
is directed almost exclusively towards smaller scales. This causes a front to form, beyond 
which no excitation has yet taken place. As this front penetrates the small scale regions, it 
leaves behind a power law distribution of the energy, whose exponent is close to a = 1.12. 
Eventually the front feels the viscosity, which happens at the Kolmogorov length rj, estimated 
from the initial conditions and the viscosity. Since the energy is now dissipated instead of 
transferred, the front stalls, and an equilibrium of inertial transfer and energy dissipation is 
established. Now there is also significant backflow of energy, and a transfer towards smaller 
scales is observed only on the average. Thus the profile gradually reduces in steepness and 
converges to the familiar Kolmogorov form |TT],^, with a scaling exponent close to the 
classical value of 2/3. Since there is no energy input, a truly stationary state cannot be 
established, and all excitations will decay to zero in the infinite time limit. This however 



will happen on much longer time scales than seen in Fig. |l[ 

From this we observe that a Kolmogorov state develops from the interplay between 
singular motion and viscous dissipation. We will now concentrate on the early time behavior, 
where viscosity is not yet important. Our aim is to explain the value of the scaling exponent 
a > 2/3, and to find the structure of the singular front. 

III. THE EULER SINGULARITY 

Here we describe the evolution of very long cascades, where the viscosity has been turned 
off. Since we only look at a single trajectory leading up to the singularity, and no statistics 
have to be accumulated, we can easily afford to simulate 100 levels, corresponding to 30 
orders of magnitude in scale. This will allow us to identify the scaling behavior of the 
singularity unambiguously. 

The result of a simulation, where again only the level £ = is excited, is shown in Fig.|^. 
Note that the axes are logarithmic, so the small scales contain only a very small fraction of 
the energy. 
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FIG. 2. The energy cascade for the REWA model with 100 levels. Energy is concentrated in 
the top level initially. The energy drops to zero at a finite level number. This position moves in 
time intervals which scale geometrically as the singularity is approached. 
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The distribution of energies Ei looks very similar to the previous figure, except that the 
absence of viscosity allows the cascading to continue indefinitely. Since the length scale 
associated with level i is r = 2~^L, the energy spectrum behind the front is a power law 

Eer^r" ,a = 1.12 ±0.01 . (5) 

At any given time, the energy drops to zero at the front. The scale where this happens 
thus represents the smallest excited scale, which goes to zero at a finite time t* which depends 
on initial conditions. Thus one expects sufficiently high derivatives of the velocity field to 
blow up as the time difference from the singularity 

T = t* -t (6) 

goes to zero. Indeed, plotting the smallest excited scale Vc as a function of r, one again finds 
a power law 
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FIG. 3. The level where the energy drops to zero, as function of the time distance from the 
singularity. The exponent /3 is related to a by the scaling relation a = 2 — 2/(3. 

r^^T^ ,/3 = 2.27 ±0.01 , (7) 

see Fig. ^ A relation between the two power laws and (|^ is established by the following 
argument: the time scale on which the singularity is moving must be r itself, so a typical 
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energy at the front is {rdrY ~ ^1^^ , Comparing this with the scahng law one finds 
the scahng relation 

2(/3-l) 



a 



/9 



(8) 



which is obeyed precisely by the values found for a and (3 numerically. 

The typical velocity is from (|^) expected to go down like r"/^ in scale. Thus the vorticity 
of a shell £, defined by 



E kxu(k,t) 



ke/c. 



(9) 



behaves like r"/^ ^ = r and reaches its maximum near the front, as seen in Fig. 




Since uj has units of inverse time, this maximum diverges like r^^, which we confirmed 
numerically. This means the singularity observed in the REWA model is consistent with the 
criterion by Beale, Kato, and Majda that the maximum of the vorticity should diverge 



at least as fast as r ^ for a true Euler singularity. Because the exponent of the vorticity is 
known, we used the scaling relation 

max{cj^} ~ {t* - ty^ (10) 
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to fit the value of tlie singular time t*. 

Next we look at the possible influence of initial conditions on the singularity. In Fig. ^ 
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FIG. 5. The singularity for two different initial conditions. If the exponent is smaller than 2 
initially, only the front moves. If it is 2, all scales move on the same time scale initially. 

we chose the energy to have a nontrivial distribution at the initial time. For this distribution 
we chose two power laws, one with an exponent smaller than a, the other larger. As seen 
in the figure, in both cases the solution settles on the same slope, with the same universal 
shape at the front. We can thus conclude that the singularity is universal except for the 
value of the singular time t* and the energy scale. This is of course only true apart from 
small fluctuations of the energy. These result from the complicated chaotic motion of Fourier 
modes underlying the excitation of small scales. However at the front one is very far from 
equilibrium, so fluctuations in the energy transfer are small compared with its absolute 
value. 

There is another interesting observation to be made in Fig. ^ which hints to the observed 
universality. If the exponent of the energy distribution is 0.8 initially, only the front of 
the distribution moves, the contribution from larger scales remains static. This is because 
Tr = {r'^/Ey^'^ = r^~"/^ represents a local time scale. Hence as long as a < 2, only the 
smallest available scale moves, since it has the shortest time scale. The limiting case a = 2 
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is the other initial distribution given in Fig. ^, and indeed it now evolves on all length scales. 
But as soon as the newly formed front overtakes the old one, it again only grows from its 
front, since the slope behind it is now smaller than 2. Hence in each case universality results 
from the fact that growth is determined only from the local properties of the front. 

Having seen that the characteristic length and time scales of the singularity behave like 
power laws, we see next whether the whole sequence of profiles can be rescaled to fall onto 
a single master curve. To that end, we first nondimensionalize length and time. If L is a 
length where the energy already has its scaling form, and Eq is the energy on that scale, one 
can introduce the nondimensional quantities 

r = r/L, Ee = E,/ E,, i = t{E^/L''f'^ . (11) 

Using these, the energy is expected to scale like 

E,(f)=f"$(0 ,e = ^ , (12) 
where $ is a universal function. In Fig. |^ 
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FIG. 6. The scaling function of the energy ^(^) determined from simulations of the REWA 
cascade using the rescaling ([l2|). It is well fitted by a functional form proposed originally for 
singularities of the full Euler equation. 
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the energies at different distances from the singularity are superimposed according to (f2|). 
The values of the energy at different levels is marked by crosses. Allowing for some fluctua- 
tions, the scaling relation is obeyed very well. Owing to the rescaling (|1T]) , (|I2[) , $ asymptotes 
to 1 as ^ CX3, and the collapse is the same for all initial conditions. Finally we note that 
$ is fitted very well by 

$(0 = exp(-5 • 10-VO , (13) 

which is a functional form of the energy spectrum proposed by Brachet et al. [@]. We will 
comment on the relation between our work and in the discussion. 



IV. TOWARDS AN ANALYTICAL DESCRIPTION 

We now develop an analytical theory for the form of the singularity which will also explain 
its universality, i.e. independence of initial conditions. As a basis we use a simple model 
for the energy transfer |T^, which has been used before to describe the stationary state of 



developed turbulence |T^. The energy transfer is split into two parts, one deterministic, the 
other stochastic: 

T,^,+i it) = T^lt 1 it) + Ti'XI it) ■ (14) 

Once T£^£+i(t) is specified, conservation of energy (^ results in an equation of motion 
for the energy. The deterministic part expresses the tendency of the cascade to establish 
equipartition of energy between its members. The stochastic part represents the chaotic 
mixing of the individual Fourier modes. If in addition we assume that both Tfi^l^^-^{t) and 
Ti^i^iit) only depend on the neighboring values of the energies, one ends up with the 
expressions 



T^tt.^t) = (Ef (t) - E%l{t)) , (15a) 

/ 9(^+1) \ -"-/^ 

Ttitht) = R{—^] iEem,^i{t)f%^^{t) . (15b) 
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The powers appearing in (|T^) are derived from dimensional considerations, and ^ represents 
a Gaussian white noise with {(.tit)) = and {C,£(t)C,i'{t')) = 26eii6(t — f). We use Ito's 
definition in equation ( |15b| ). Together with (|T5|) is a Langevin equation for the motion 
of a cascade, so we will refer to it as the Langevin model. It has been shown that the model 
given by (^,(|TB|) exhibits multifractal scaling in a stationary turbulent state, and anomalous 



scaling exponents can be calculated analytically [|l6l. At the same time it gives an excellent 
description of the turbulent state of the REWA cascade [|14 |. 
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FIG. 7. The same scaling function $ as in Figure |6|, but for the Langevin model ([^,([151). The 
crosses represent the full model, the full line corresponds to the deterministic case R = 0. 



Figure ^ shows a simulation of the model equations (^, ([151) at zero viscosity, with energy 
concentrated in the largest scale initially. The profiles have been rescaled as in (|12D, which 
again leads to a collapse very similar to that of Fig. The exponent a = 1.24 is slightly 
larger than that found for the REWA cascade. The values of the free parameters D and R 
in ([T5|) are taken from [Q: 



L) = 6.5-10~2 ,i? = 4.4-10~2 . (16) 

The amplitude D is of particular significance, since it determines the effectiveness of the 
energy transfer. If D gets larger, the front reaches a given scale at an earlier time, as we are 
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going to see in more detail below. 

On the other hand, the noise strength R is insignificant for the formation of the singu- 
larity, as fluctuations are small, in agreement with the result found for the REWA cascade. 
This is to be expected since the motion of the singularity is dictated by the front which is 
very steep. Thus the deterministic part ( |15a| ), which consists of an energy difference will 



dominate the stochastic part (|15bD . We verified this by putting R to zero, leaving everything 
else unchanged. The result is shown as the solid line in Fig. |^, which is a perfect fit to the 
fluctuating data of the stochastic cascade. 

Our first approximation will thus be to include only the deterministic part of the Langevin 
model in our analytical description. A second approximation is of a more technical nature. 
It is seen in Fig. |^ that the sequence of level energies form a reasonable approximation 
of a continuous curve. This motivates us to pass to a continuum limit, where the ratio of 
length scales between two levels approaches one, leading to a less cumbersome description 
in terms of differential equations. We introduce A as the ratio of length scales, which means 
that the length scale on level i is r = X~^L. As A approaches 1, Ei can be written as a 
continuous variable E = E{r). Replacing 2~^L by A~^L in and and performing the 
limit A — s> 1, one ends up with 

dtE{r, t) = DrdlE^/^ (r, t) , (17) 

where D = Z}(lnA)^ is a rescaled coupling constant. The equation of motion ([T7|) is the one 
our subsequent analytical description is based on. 

We now look for self-similar solutions of ([17D of the form 

^(f,f) = f"$(e) ,e = r/f^ , (18) 

in direct analogy to (|12|) . Plugging this into the equations of motion, we find that the explicit 
dependence on r is eliminated by demanding that /3 = 2/(2 — a), so we recover the scaling 
relation (^. As a result, we are left with a similarity equation which depends on ^ = r/f^ 
alone: 
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I — a III z 

+e^''\m\d ■ (19) 



In the following section we will show that (|T^) possesses a unique physical solution, which 
fixes both the similarity function $ and the exponent a. 



V. SELECTION 

Since (|19D is of second order, one needs two initial conditions to uniquely fix the solution. 
As noted earlier, $ asymptotes to a constant at infinity. Since scales have been normalized 
according to (|11|), this constant is one, leaving us with the boundary conditions 

$'(0 ^ 

The coupling strength D can be eliminated by the transformation 

^^^2/(a-2)^ . (21) 

This means that for large coupling strengths the position of the front moves towards larger ^. 
Thus a given length scale is reached earlier, as to be expected on physical grounds. For our 
discussion of universal solutions we will consider the equation in the independent variable 
where D has been eliminated. 

To find explicit solutions, we expand the similarity equation around C = oo. Using the 
boundary conditions, this leads to an asymptotic expansion 

oo 

*(C) = E^^C-^^'' ' ^0 = 1 , (22) 
i=n 

where the coefficients hi depend on a alone. Given a sufficiently large Qnit, ( p2[ ) can be used 
to generate an initial condition at Cmit- This initial condition then allows to integrate the 
similarity equation numerically towards small C,. Hence there is a unique solution for each 
a, while we expect the partial differential equation (0) to select a unique a. To understand 
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oo 



(20) 



this, we now look at the behavior of solutions for different a. Three cases arise, which are 
shown in Fig. ||. 
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FIG. 8. Solutions of the similarity equation (19) for three different a. For a greater than a 
critical value Oc = 1.12045, the tip of the front becomes fractal. The inset shows a blow-up of this 
behavior on a logarithmic scale. 



If a is smaller than a critical value ac = 1.12045 ± 10~^ the profile ends in a sharp 
front, similar to the front observed in simulations both of the REWA and the Langevin 
cascade. The profile is zero below a front position (q, for ( > (q ^ behaves asymptotically 
like $ ~ (C — Co)^''^- This corresponds to a local expansion of the form 

oo 

m = T.^^i(-(of^'^^' ■ (23) 

If on the other hand a is larger than ac, the front becomes unstable and levels off to 
form a flat plateau. As shown in more detail in the appendix (see also the inset in Fig. 

this plateau dips down again at a smaller value of (, only to form another plateau as 
the second front becomes unstable. This process repeats itself, to form a fractal tip which 
asymptotes to zero. Clearly this is not an acceptable solution, at the very least because it 
would correspond to energy being transported instantly across all levels. 

This indicates that there is something special about solutions at a = ac, which separates 
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the region of sharp and fractal fronts. Indeed, a more careful analysis, which is detailed in 
the appendix, reveals that at the critical a the asymptotics at the front position is now 

oo 

$(C)=E«^(C-Cor^ • (24) 

But although this amounts only to a slight difference in the appearance of the fronts in Fig. 
^, the solution at Uc is the one which is selected. This comes from an argument similar to 
that advanced for front propagation into an unstable medium |^8|. Indeed, on a logarithmic 



scale, i.e. by level number, the self-similar solution ([TsD corresponds to a front propagating 



at a constant speed a. In our problem, the situation is actually reverse to that of |2^: of 
all possible solutions, the one with the highest speed will eventually take over, while the 
slow solutions are left behind. On the time scale set by the front, they no longer move, and 
thus drop out of the problem. This explains the universality observed earlier: independent 
of initial conditions, only one front with a given exponent a is observed. We also checked 
the validity of our selection argument directly, by simulating the Langevin cascade (^), (|15|) 
with i? = for smaller and smaller scaling factors A. Extrapolating to A = 1, we were able 
to confirm the value of Uc to five decimal places. 

In Table | we summarize some of the values of the exponent a obtained for different 
cascades. The REWA cascade with 26 modes, which we considered throughout this paper, 
is called "small cascade" here, to distinguish it from another mode selection with 74 modes. 
The Langevin cascade with the same scale factor A = 2 as the REWA cascades gives a 
somewhat larger value. However, the overall variation of the exponent a is only in the 
order of 10%. This strongly supports the claim that our analytical theory has captured the 
physical mechanism behind the selection of a singular front for all the models considered. 

Another important quantity is the position of the front, which for the similarity solution 
is found to be 

f, = (25) 

Once again we see that for large D a given length is excited at earlier times. The coupling 
D is the only parameter to be determined for a comparison between theory and simulation 



a 



small cascade, N=26 
large cascade, N=74 
Langevin cascade, A = 2 
similarity equation 



1.12 ±0.01 
1.17 ±0.01 
1.24 ±0.01 
1.12045 ± 10"^ 



TABLE L Compilation of different values of the exponent of the shell energies a. The first two 
values are from simulations of the REWA cascade with two different wave vector sets. The third 
line refers to the model equations (^) , (^) , the last line is the result of our similarity theory. 

of the REWA cascade. We adjusted D such that the average energy in a turbulent state 



agrees with the value determined for the REWA cascade flj]. Thus we are able to predict 
$ without adjustable parameters. The comparison between the solution of the similarity 
equation at a = ac and the numerical simulation of the REWA cascade is shown in Fig. 0. 
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FIG. 9. Comparison between the similarity theory and and the numerical simulation of the 
REWA cascade. The free parameter D of the similarity equation has been determined from an 
independent measurement of the turbulent state. 



The agreement is quite good, although there is some difference in the shape of the front. 
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It makes sense that for this shape, which describes the excitation of Fourier modes, the 
detailed couphng structure of modes matters to some degree. This shows that the formation 
of a singularity in a very complicated system of coupled Fourier modes is indeed described 
by a very simple physical principle: the tendency of the cascade to establish equipartition 
of energy. 

VI. DISCUSSION 

Roughly speaking, two different methods have been used to numerically treat the for- 
mation of singularities in the Euler equation. One relies on a Fourier representation of the 
velocity field 0,^] , the other tries to track the vorticity in real space . Given a simulation 
with perfect resolution, the results should be the same, but in effect there are two different 
physical pictures underlying the two approaches. A spectral approach is preferable if a sin- 
gularity generically results from the interaction of Fourier modes, as suggested originally by 
Onsager A spatial picture would show a uniformly (multi)-fractal structure. 

Vorticity dynamics, on the other hand, would be more useful if a singularity results from 
a specific spatial structure, like the meeting of two antiparallel vortex tubes P]. Other, more 
complicated structures have been proposed in two-dimensional, axisymmetric flow 

Since the spatial resolution of the REWA cascades studied here is quite low, it is hard 
to speculate about the physics of vortex tubes. We must rely on the spectral picture being 
relevant, which was analyzed in detail by Brachet et al. 0. They find that the large 
wavenumber part of the energy spectrum is well described by 

E(A;,t) = aA;-"We-2^W'= . (26) 

This is very similar to the form of the energy distribution observed for our case. However, an 
important difference is that (^) predicts excitation of all scales, while in the REWA model 
the energy drops to zero. This difference, which becomes noticeable for energies smaller 
than those shown in Fig. ^, comes from the absence of non-local interactions in the REWA 
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model, which would excite small scales instantly. Since the energy Eg of a shell represents 
the energy spectrum integrated over an octave in wavenumber, n{t) is to be identified with 
the exponent a + 1. If in addition 6{t) was chosen to behave like t^, one would end up with 
our scaling form ([T8|). 

However, although n{t) converges to a value close to 4 0, this is much larger than the 
value a+1 = 2.12 we find. Other workers find n close to 3 asymptotically. This is 
closer to our value, but still at the boundary where the largest scales would move fastest, 
according to our estimate Tr = r^^"/^. Our simulations are self-consistent since for a > 2 
nonlocal interactions would become dominant which have not been included in 



our description. It would of course be of great importance to investigate whether local or 
non-local interactions are dominant for the full Euler equation. 

A second significant difference between the results of and our work is that 6{t) does 
not go to zero in finite time, but rather behaves like 

S{t) = Soe-'/^ . (27) 

This corresponds to a singularity only in infinite time. All results of were later confirmed 
in 1^ using greater resolution and more general initial conditions. But of course there is also 
a good possibility that if greater resolution allows to continue the simulation still further, a 
crossover from (^71) to a finite time singularity is seen. 

A definite answer whether singularities of the three-dimensional Euler equation exist 
can only be given by considering the full equations. Nevertheless, our study shows that 
structures exist in inviscid fiow which show blow-up in finite time. They are consistent with 
the divergence of the vorticity like r~^, and the existence of a local cascade. The scaling 
structure of the singularity found here can thus serve as a guideline for further studies of the 
full equations. The more general structures possible in the fully resolved fiow can of course 
be more singular than ours, so we expect the structures appearing in our mode-reduced 
systems to be dominated by even more singular contributions. Conversely, there is also the 
possibility that the modes not taken into account in our study will interact with the reduced 
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mode system to keep it from becoming singular. One possible way to study this would be 
to selectively take non-local interactions into account, to find out what modifications they 
imply for our analysis. 

A particularly intriguing aspect of the present work is the novel way the scaling exponent 
is selected as the "speed" of a marginally stable solution. This adds another variety to 
the existing mechanisms which determine the scaling exponents of singularities which are 
not determined from dimensional arguments. Other selection mechanisms for this scaling 
behavior of the second kind (in Barenblatt's terminology) are found in and in 



Formally, the selection mechanism is quite similar to that of marginally stable solutions 



of equations of the form [28 



dt<j) = dl<p + F{<f)) . (28) 

However, solutions to this equation do not drop to zero at a finite value of x but rather 
decay exponentially away from the front. For this reason we are not able to repeat the 
linear stability analysis presented for example in since the front of our solution is very 



steep. In that respect it is more similar to solutions of the porous medium equation 
which in one dimension reads 

dtu = dlu"^ . (29) 

Solutions to this equation form a front which drops to zero, like ours. However, although 
for m = 3/2 ([29| ) looks quite similar to (p!?!), the factor r in (|1^ represents a very serious 



complication. This is because its value at the front goes to zero at the singularity, giving a 
very singular diffusion constant. Thus again we are not able to carry over the mathematically 
rigorous results known for the porous medium equation. Still we believe that the combination 
of numerical evidence and analysis of the similarity equation conclusively demonstrates that 
the marginal solution is indeed the one selected for our equation ([T7|). This adds to the 
generality of the marginal stability concept. 

Finally, we would like to stress the fact that the cascades studied in this paper give a 
coherent picture of both Euler singularities and fully developed turbulence. In fact the same 
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is true for a scalar shell model recently studied by Dombre and Gilson ||3^ . They also found 



unique singular solutions of the inviscid equations, but whose spectrum is less steep than 
Kolmogorov's. In addition, they propose a connection between inviscid singularities and 
intermittent fluctuations. We started to explore this connection in the second section, but 
this has to be pursued further. We suspect the situation will be quite different from that 
of the scalar shell model, because our singular spectrum is steeper than that of a stationary 
cascade. In particular, it would be interesting to study the interaction between the most 
singular solutions and viscosity, leading to a turbulent state. Since the most singular solution 
with a = etc is regularized by viscosity, other solutions with a < ac will become relevant as 
well. So it may be the whole spectrum of singular solutions which is relevant to the turbulent 
state. 
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APPENDIX: SIMILARITY SOLUTIONS 

Here we discuss the transition of solutions of the similarity equation (|l^) from regular 
to fractal tips in more detail. The similarity variable is rescaled according to C = -D^''""^^. 

The invariance of ([T9|) under scale transformations ( fi(, $ fi'^~'^^ can be used to 
transform it to a first order equation 



v'(u) 



u^/'^v{u) 
with 



3 6-3a ^ ^ ^ ' 2 ^ \ 



(Al) 



v{u) = Cu'{0 . (A2) 
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Primes always refer to derivatives with respect to the argument. It follows from ( |A2|) that 
u goes to zero both for C going to infinity and $ going to zero. Thus all solutions start out 



at M = 0, corresponding to ( = oo, shown by the full line in Fig. |10|. To form a regular tip, 
they must return to u = (dashed and dot-dashed lines), which is the case for a < ac- For 
a > ac, on the other hand, v{u) enters a limit cycle, which corresponds to the fractal tip 
(dotted line). 
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FIG. 10. Solutions of the first-order version of the similarity equation. The three different types 
of solutions correspond directly to Fig. |^. The fractal tip, observed for a > ac corresponds to a 
limit cycle. 



Equation (|A1|) becomes singular for u = or f = 0, thus only local solutions exist either 



in the domain < u < Umnr or < u < Umnr with viu. 



^max/min i 



0. To form a complete 



solution, local solutions have to be matched at the points indicated in Fig. |T^. 
Three different expansions exist around the point = 0: 

oo 

7=1-11 
' ' 2'2 ' 



(A3) 



which corresponds directly to the three expansions (^), (^31) , and (p^ of $. The coefficients 
cf'\ which describe the behavior of $ at infinity, are determined recursively from 
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4 ^45, + 5c.^i\ + E(i + |jcy^4Vij • (A4) 

The case 7 = 1/2, which corresponds to a tip of the form $(C) = Ccr°(co^^^'')^(C ~ Co)^/4 is 
determined by 

Co 



3(2 -a 

(1/2) _ -26a + 65.2 + 154^1^72 + (3/2 + 3j7 4) '^cjif 
2/(2 -a) -3 (1 + ^/4) ci'/') 



(A5) 



This is the critical case. On the other hand, solutions with a < ac have a tip $(C) = 
Co~"(3co ^^^''/2)*'^^^''(C~Co)^^'^- In that case the constant Cq is open and the higher order 
coefficients are calculated from 



J-l/2) _ ^ , 2 (_i/2) ,.(-1/2),^ Y^/Q./^N J-1/2) f-1/2) 



25.5 - 65.6 + ^ diZ^'^ - 154:3''V2 - E(3j74)c} 



c 



o- (-1/2) P'"*^ ""^t) ' 9 -""--1-3 A^K'^J/^J^j '-t-j 



(A6) 



All expansions ( [A4| )-([A^) are asymptotic in nature, but are extremely good everywhere 



except in a small neighborhood of Umax- To do the matching at Umax, equation ( [Al[ ) has to 
be integrated numerically. 

The local behavior at Umax is 

V{U) = C {Umax - uf' with C = ± (^8{Umax " ^^I'j) , (AT) 

where the solid line of Fig. |10] corresponds to negative C, the other branches to positive 
C. At a given a, for each of the three branches 7 = 1, —1/2, and 1/2 one can extrapolate 
to the asymptotic behavior ( [A7| ) to determine Umax- A matching of the branch 7 = 1 on 
the branch 7 = 1/2 is only possible for one a, which is the critical a = ac = 1.12045. For 
7 = — 1/2 the free parameter Cq~^^'^^ can be used to match Umax, and this turns out to be 
possible only for a < ac- 

If a > ac, there is no branch which returns to m = 0, and the solution has to be 
continued numerically through the matching point. The resulting branch goes to f = at a 
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finite u = Umin, where it can be matched onto the original solution with 7 = 1. This means 
the solution ends up on a limit cycle. From the definition of u (cf. ( |A2| )) it is clear that the 
sequence of tips of $(C) approximates a power law $ ~ C^~") &s shown in the inset of Fig. 
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